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November 2, 2023 at 22:24 


1. Intro. This program is part of a series of “exact cover solvers” that ’m putting together for my own 
education as I prepare to write Section 7.2.2.1 of The Art of Computer Programming. My intent is to have 
a variety of compatible programs on which I can run experiments, in order to learn how different approaches 
work in practice. 

Instead of actually solving an exact cover problem, DLX3-PRE is a preprocessor: It converts the problem 
on stdin to an equivalent problem on stdout, removing any options or items that it finds to be unnecessary. 

The program DLX-PREallows preprocessing of XCC problems. This program adds support for MCC 
problems. The extensions were done by Filip Stappers. 

The input format for this program is the one described in DLX3. Most of the code below, like the 
description above, has been cribbed from DLX3 and DLX-PRE, with minor changes. The output is similar 
to the format described in DLX3-PRE, but will now also output the multiplicities of the items. 

After this program does its work, it reports its running time in “mems”; one “mem” essentially means a 
memory access to a 64-bit word. (The given totals don’t include the time or space needed to parse the input 
or to format the output.) 

Here is the overall structure: 


7##define o memst+ /* count one mem */ 

7##define oo mems += 2 /* count two mems */ 

#define 000 mems += 3 /* count three mems */ 

#define O "%" /* used for percent signs in format strings */ 

#tdefine mod % /* used for percent signs denoting remainder in C */ 

7##define maz_level 500 /* at most this many options in a solution */ 

#define maz_cols 10000 /* at most this many items */ 

7##define maz_nodes 10000000 /* at most this many nonzero elements in the matrix */ 


#tdefine bufsize (9 * maz_cols + 3) /* a buffer big enough to hold all item names «/ 


#include <stdio.h> 
#include <stdlib.h> 
#include <string.h> 
#include <ctype.h> 
#include "gb_flip.h" 
typedef unsigned int uint; /* a convenient abbreviation */ 
typedef unsigned long long ullng; /* ditto */ 
( Type definitions 4); 
(Global variables 2); 
(Subroutines 8); 


int main(int argc,char *argv[]) 
{ 
register int c, cc, dd, i, 7, k, p, pp, q, 99, 7, rr, Trr, t, wu, x, cur_node, best_itm, stage; 
(Process the command line 3); 
(Input the item names 13); 
(Input the options 1s ); 
(Tighten bounds of the primary items 22); 
if (vbose & show_basics) (Report the successful completion of the input phase 23); 
if (vbose & show-_tots) {Report the item totals 24); 
imems = mems,mems = 0; 
(Reduce the problem 29); 
finish: (Output the reduced problem 43); 
done: if (vbose & show_tots) {Report the item totals 24); 
all_done: if (uvbose & show_basics) { 


} 


2. 


INTRO DLX3-PRE 
forintf (stderr, 
"Removed,,"O"doption"O"s,,"O"d item"O"s and tightened,"O"dbound"O"s", 
options_out, options_.out =1?"":"s", itmsout, itms.out =1?"":"s", bounds_tightened , 
bounds_tightened = 1? "":"s"); 
forintf (stderr, "after,"O"1llut+"O"1llu_mems,", imems, mems); 
forintf (stderr, "\"O"dyround"O"s.\n", rnd, rnd =1?"":"s"); 


} 


options on the command line: 


1 


You can control the amount of output, as well as certain properties of the algorithm, by specifying 


e ‘v(integer )’ enables or disables various kinds of verbose output on stderr, given by binary codes such as 
show_choices; 
e ‘d( integer )’ to sets delta, which causes periodic state reports on stderr after the algorithm has performed 
approximately delta mems since the previous report (default 10000000000); 
e ‘t (positive integer )’ to specify the maximum number of rounds of option elimination that will be at- 
tempted. 


e ‘T( integer ) 


? 


clause, but doesn’t ruin the integrity of the output). 


7##define show_basics 1 /* vubose code for basic stats; this is the default */ 
#define show_choices 2 /* vbose code for general logging */ 


7##define show_details 4 /* vubose code for further commentary */ 
#tdefine show_orig_nos 8 /* vbose code to identify sources of output options */ 
7##define show_tots 512 /* vubose code for reporting item totals at start and end */ 
#define show_warnings 1024 /* vbose code for reporting options without primaries */ 
(Global variables 2) = 

int vbose = show_basics + show_warnings; /* level of verbosity */ 

char buf [bufsize]; /* input buffer «/ 

ullng options; /* options seen so far */ 

ullng imems, mems; /* mem counts */ 


ullng thresh = 1000000000; /* report when mems exceeds this, if delta £0 */ 
ullng delta = 10000000000; /* report every delta or so mems */ 

ullng timeout = *1ffffffff{fffffff; /* give up after this many mems */ 

int rounds = maz_nodes; /* maximum number of rounds attempted */ 

int options_out, itms_out, bounds_tightened; /* this many reductions made so far */ 


See also sections 6 and 30. 


This code is used in section 1. 


sets timeout (which causes abrupt termination if mems > timeout at the beginning of a 
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3. If an option appears more than once on the command line, the first appearance takes precedence. 


(Process the command line 3) = 
for (j = argc —1,k =0; 7; j——-) 

switch (argu [j][0]) { 
case ’v’: k |= (sscanf (argu[|j] + 1,""O"a", &ubose) — 1); break; 
case ’d’: k |= (sscanf (argu[j] + 1,""O"11d", & delta) — 1), thresh = delta; break; 
case ’t’: k |= (sscanf (argv[j] + 1,""O"a", &rounds) — 1); break; 
case ’T’: k |= (sscanf (argv|j] + 1,""O"1lda", &timeout) — 1); break; 
default: k = 1; /* unrecognized command-line option */ 


} 

if (k) { 
fprintf (stderr, "Usage: ,"O"s,[v<n>] |, [d<n>] ,, [t<n>] |, [T<n>] .<foo.d1x,>,bar.d1x\n", argu (0]); 
exit (—1); 

} 


This code is used in section 1. 
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4. Data structures. Each item of the input matrix is represented by a item struct, and each option is 
represented as a list of node structs. There’s one node for each nonzero entry in the matrix. 

More precisely, the nodes of individual options appear sequentially, with “spacer” nodes between them. 
The nodes are also linked circularly within each item, in doubly linked lists. The item lists each include a 
header node, but the option lists do not. Item header nodes are aligned with a item struct, which contains 
further info about the item. 

Each node contains four important fields. Two are the pointers up and down of doubly linked lists, already 
mentioned. A third points directly to the item containing the node. And the last specifies a color, or zero if 
no color is specified. 

A “pointer” is an array index, not a C reference (because the latter would occupy 64 bits and waste cache 
space). The cl array is for item structs, and the nd array is for nodes. I assume that both of those arrays 
are small enough to be allocated statically. (Modifications of this program could do dynamic allocation if 
needed.) The header node corresponding to cl{c] is nd[c]. 

Notice that each node occupies two octabytes. We count one mem for a simultaneous access to the up 
and down fields, or for a simultaneous access to the itm and color fields. 

This program doesn’t change the itm fields after they’ve first been set up, except temporarily. But the 
up and down fields will be changed frequently, although preserving relative order. 

Exception: In the node nd|[c] that is the header for the list of item c, we use the itm field to hold the 
length of that list (excluding the header node itself). We also might use its color field for special purposes. 
The alternative names len for itm and aux for color are used in the code so that this nonstandard semantics 
will be more clear. 

A spacer node has itm < 0. Its up field points to the start of the preceding option; its down field points 
to the end of the following option. Thus it’s easy to traverse an option circularly, in either direction. 

Spacer nodes are also used within an option, if that option has been shortened. The up and down fields 
in such spacers simply point to the next and previous elements. (We could optimize this by collapsing links, 
for example when several spacers are consecutive. But the present program doesn’t do that.) 


#define len itm /* item list length (used in header nodes only) */ 
#define aux color /* an auxiliary quantity (used in header nodes only) «/ 
(Type definitions 4) = 

typedef struct node_struct { 


int up, down; /* predecessor and successor in item */ 

int itm; /* the item containing this node */ 

int color; /* the color specified by this node, if any */ 
} node; 


See also section 5. 


This code is used in section 1. 


5. Each item struct contains five fields: The name is the user-specified identifier; next and prev point to 
adjacent items, when this item is part of a doubly linked list. bound and ubound are the lower and upper 
bounds for the multiplicity of this item. The lbound and ubound fields of secondary items are not used. 
We count one mem for a simultaneous access to the prev and nezt fields, or for a simultaneous access to 
Ibound and ubound. 
(Type definitions 4) += 
typedef struct itm _struct { 


char name[8}; /* symbolic identification of the item, for printing */ 
int prev, neat; /* neighbors of this item */ 
int lbound, ubound; /* allowed multiplicities */ 


} item; 
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6. (Global variables 2) += 
node nd[maz_nodes]; /* the master list of nodes */ 
int last_node; /* the first node in nd that’s not yet used «/ 
item cl[maz_cols + 2]; /* the master list of items «/ 
int second = maz_cols; /* boundary between primary and secondary items */ 
int last_itm; /* the first item in cl that’s not yet used */ 


7. One item struct is called the root. It serves as the head of the list of items that need to be covered, 
and is identifiable by the fact that its name is empty. 


##define root 0 /* cl[root] is the gateway to the unsettled items */ 


8. An option is identified not by name but by the names of the items it contains. Here is a routine that 
prints an option, given a pointer to any of its nodes. It also prints the position of the option in its item. 

This procedure differs slightly from its counterpart in DLX2: It uses ‘while’ where DLX2 had ‘if’. The 
reason is that DLX-PRE sometimes deletes nodes, replacing them by spacers. 


(Subroutines 8) = 
void print_option(int p, FILE «stream) 
{ 
register int k, q; 
if (p < last_itm V p > last_node V nd[p].itm <0) { 
forintf (stderr, "Illegaloption,"O"d!\n", p); 
return; 


for (q=p;; ) { 
fprintf (stream, "\"O" .88", cl[nd|q].itm].name); 
if (nd[q|.color) fprintf (stream,":"O"c", nd[q|.color > 0? nd[q].color : nd[nd[q].itm].color); 
qn; 
while (nd|q].itm <0) ¢ = nd{q].up; 
if (¢ =p) break; 
} 
for (q = nd[nd[p].ttm].down,k =1; q#p; k++) { 
if (q = nd[p].itm) { 
forintf (stream, "(?)\n"); return; /* option not in its item! */ 
} else gq = nd[q].down; 
f 
forintf (stream, "WC"O"dyofy"O"d) \n", k, nd[nd[p].itm].len); 
} 
void prow (int p) 
{ 
print_option (p, stderr ); 
} 
See also sections 9, 10, 11, 27, and 28. 


This code is used in section 1. 
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9. Another routine to print options is used for diagnostics. It returns the original number of the option, 
and displays the not-yet-deleted items in their original order. That original number (or rather its negative) 
appears in the spacer at the right of the option. 
(Subroutines 8) += 

int dpoption(int p, FILE «stream) 


register int q, c, first; 
first = 1; 
for (p—-; nd[p].ctm >0V nd[p].down < p; p—) ; 
for (¢=p+1; ; q+) { 
c= nd{q].itm; 
if (c <0) return —c; 


if (c>0) { 
if (first) fprintf (stream,"\"); 
else first = 0; 


fprintf (stream, ""O" .88", cl[c].name); 
if (nd|q|.color) fprintf (stream,":"O"c", nd[q].color); 
} 


} 
} 


10. When I’m debugging, I might want to look at one of the current item lists. 


(Subroutines 8) += 
void print_itm (int c) 
{ 
register int p; 
if (c < root Vc > last_itm) { 
forintf (stderr, "Illegaljitem,"O"d!\n",c); 
return; 


fprintf (stderr, "Item,"O".8s", cl[c].name); 
if (c < second) { 

if (cl[{c|.lbound & cl[c|.ubound V cl{c].ubound F 1) 

fprintf (stderr, "{C"'O"d,"O"da)", cl[c].lbound, cl[c].ubound); 
forintf (stderr ,", Jlength,,"O"d, yneighbors,."O".8syjand,"O".8s:\n", nd[c].len, 
cl[cl[c].prev].name, cl{cl{c].next].name); 

} else fprintf(stderr," , length,"O"d:\n", nd[c].len); 
for (p = nd[c].down; p> lastitm; p= nd[p|.down) prow(p); 
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11. Speaking of debugging, here’s a routine to check if redundant parts of our data structure have gone 
awry. 
7##define sanity_checking 0 /* set this to 1 if you suspect a bug */ 
(Subroutines 8) += 
void sanity (void) 


{ 


register int k, p, q, pp, 9, ¢; 
for (q = root, p = cl[q].neat; ; q=p,p = cl[p|.next) { 
if (cl[p|.prev # q) fprintf(stderr, "Badyprev.fieldjat itm ,"O".8s!\n", cl[p].name); 
if (p = root) break; 
(Check item p 12); 
} 
} 


12. (Check item p 12) = 
for (qq = p, pp = nd{qq].down,k = 0; ; aq = pp, pp = nd[pp].down, k++) { 
if (nd[pp].up 4 qq) fprintf (stderr, "Bad upfield at modey"O"d!\n", pp); 
if (pp = p) break; 
if (nd[pp].itm 4 p) fprintf(stderr, "Bad itmfieldjat_node,"O"d!\n", pp); 


f (nd[p|.len #k) fprintf(stderr, "Badlen, field in item,"O".8s!\n", cl[p].name); 


This code is used in section 11. 
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13. Inputting the matrix. Brute force is the rule in this part of the code, whose goal is to parse and 
store the input data and to check its validity. 


#define panic(m) 


{ fprintf(stderr, ""O"s!\n"O"d:4"O".99s\n",m, p, buf); exit (—666); } 


(Input the item names 13) = 


if (maz_nodes < 2 * maz_cols) { 
fprintf (stderr, "Recompile_me:max_nodes,mustjexceed, twice max_cols!\n"); 
exit (—999); 
} /* every item will want a header node and at least one other node »/ 
while (1) { 
if (fgets (buf, bufsize, stdin)) break; 
if (0, buf [p = strlen(buf) — 1] A ’?\n’) panic("Input,lineway,jtoo,long"); 
for (p = 0; 0, isspace (buf [p]); p++) ; 
if (buf [p] =’ |? V buf [p]) continue; /* bypass comment or blank line «/ 
last.itm = 1; 
break; 


if (slast_itm) panic("No,items"); 
for (; 0, buf[p]; ) { 
(Scan an item name, possibly prefixed by bounds 14); 
(Initialize last_itm to a new item with an empty list 17); 
for (p += j +1; 0, isspace(buf |p]); p++) ; 
if (buf [pl =? 1?) { 
if (second 4 maz_cols) panic("Itemname,line,contains,,| twice"); 
second = last_itm; 
for (p++; o, isspace (buf [p]); p++) ; 
} 
} 


if (second = maz_cols) second = last_itm; 

o, cl[root].prev = second — 1; /* cl[second — 1].next = root since root =0 x/ 
last_node = last_itm; /* reserve all the header nodes and the first spacer */ 
o, nd |last_node].itm = 0; 


This code is used in section 1. 
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14. (Scan an item name, possibly prefixed by bounds 14) = 
if (second = maz-_cols) stage = 0; else stage = 2; 
startname: for (j =0; 7 < 8A (0, 7isspace (buf [p+ j])); j++) { 
if (buffp+)=?2°) { 
if (stage) panic("Illegal,‘:’ yin item name"); 
(Convert the prefix to an integer, q 15); 
r=4q, stage = 1; 
goto start_name; 
} else if (buf[p+j] = 71’) { 
if (stage > 1) panic("Illegal,‘ |? inyitem name"); 
(Convert the prefix to an integer, q 15); 
if (¢q=0) panic("Upper_bound,is,zero" ); 
if (stage =0) r=q; 
else if (r > q) panic("Lower bound exceeds upper,bound"); 
stage = 2; 
goto start_name; 


} 


o, cl[last_itm].name|j] = buf [p + J]; 


switch (stage) { 

case 1: panic("Lower,bound_without_upper,bound"); 

case 0: g=r=1; 

case 2: break; 

} 

if (7 =0) panic("Item mame,jempty"); 

if (7 = 8A 7isspace (buf [p+ j])) panic("Item name,jtoo,long"); 
(Check for duplicate item name 16); 


This code is used in section 13. 


15. (Convert the prefix to an integer, g 15) = 
for (q = 0, pp =p; pp <p+Jj; pp++) { 
if (buf [pp] < °0’ V buf [pp] > °9’) panic("Illegal,digit,,in, bound spec"); 
q=10*q+ buf [pp] — 70’; 


p= pp +1; 
while (j) cl{last_itm].name|——j] = 0; 


This code is used in section 14. 


16. (Check for duplicate item name 16) = 
for (k = 1; 0, strncmp(cl[k].name, cl[last_itm].name,8); k++) ; 
if (k < lastitm) panic("Duplicate,item,mame" ); 


This code is used in section 14. 


17. (Initialize last_itm to a new item with an empty list 17) = 
if (last_itm > maz_cols) panic("Too many,items"); 
if (second = maz_cols) 00, cl{last_itm — 1]|.neat = last_itm, cl[last_itm].prev = last_itm — 1,0, 
cl{last_itm].ubound = q, cl[last_itm].lbound = r; 
else 0, cl[last_itm].nezt = cl[last_itm].prev = lastitm; 
o, nd[last_itm].up = nd|lastitm].down = last_itm; /*x nd[last_itm].len =0 */ 
last_itm ++; 


This code is used in section 13. 
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18. In DLX1 and its descendants, I put the option number into the spacer that follows it, but only because 
I thought it might be a possible debugging aid. Now, in DLX-PRE, I’m glad I did, because we need this 
number when the user wants to relate the simplified output to the original unsimplified options. 
(Input the options 18) = 
while (1) { 
if (fgets (buf , bufsize, stdin)) break; 
if (0, buf [p = strlen(buf) — 1] A ?\n’) panic("Option line,too,,long" ); 
for (p =0; 0, isspace (buf [p]); p++) ; 
if (buf [p] =’ |? V buf [p]) continue; /* bypass comment or blank line */ 


i = last_node; /* remember the spacer at the left of this option */ 
for (pp = 0; buf |p]; ) { 
for (j = 0; j < 8A (0, >isspace (buf [p + j])) A buf |p + j] #?:?5 J++) 


o, cl[last_itm].name |j] = buf [p + J]; 
if (47) panic("Empty item name"); 
if (7 = 8A 7isspace (buf [p + j]) A buf [p+ 7] 4 ?:’) panic("Item mame, too,,long"); 
if (j < 8) 0, cl[last_itm].name[j] = ’?\0’; 
(Create a node for the item named in buf [p] 19); 
if (buf [p+ j] #’:’) 0, nd[last_node].color = 0; 
else if (k > second) { 
if ((0, isspace (buf [p + 7 + 1])) V (0, aisspace (buf [p + 7 + 2]))) 
panic("Color,must bea, single,character"); 
o, nd[last_node].color = (unsigned char) buf |p + j + 1]; 
pt+=2; 
} else panic("Primary,item must jbe,uncolored"); 
for (p += j +1; 0, isspace(buf [p]); p++) ; 


i 
if (spp) { 
if (vbose & show_warnings) fprintf(stderr, "Option,ignored,,(no_primary,items) :,,"O"s", buf); 
while (last_node > i) { 
(Remove last_node from its item 21); 


last_node —; 
} else { 
o, nd[i].down = last_node; 
last_node +; /* create the next spacer */ 
if (last_node = maa_nodes) panic("Too,many,modes"); 
options ++; 
o, nd[last_node].up =i+4 1; 
o, nd[last_node].itm = — options; 


} 
: 


This code is used in section 1. 
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19. (Create a node for the item named in buf [p] 19) = 
for (k = 0; 0, strncmp(cl[k].name, cl[last_itm].name,8); k++) ; 
if (k = lastitm) panic("Unknown,,item mame"); 
if (0, nd[k].auz > 7%) panic("Duplicate item name in this,option"); 
last.node ++; 
if (last_node = maa_nodes) panic("Too many modes"); 
o, nd[last_node|.itm = k; 
if (k < second) pp = 1; 
o,t = nd[k].len + 1; 
(Insert node last_node into the list for item k 20); 


This code is used in section 18. 


20. Insertion of a new node is simple. We store the position of the new node into nd|k].aux, so that the 
test for duplicate items above will be correct. 


(Insert node last_node into the list for item k 20) = 
o, nd[k].len = t; /* store the new length of the list */ 
nd|k].aux = last_node; /* no mem charge for aux after len */ 
o,r = ndl[k].up; /* the “bottom” node of the item list «/ 
000, nd[r].down = nd|k].up = last_node, nd[last_node].up = r, nd[last_node].down = k; 


This code is used in section 19. 


21. (Remove last_node from its item 21) = 
0, k = nd[last_node].itm; 
oo, nd|k].len-—, nd|k].auz =i—1; 
0,q = nd[last_node].up,r = nd{last_node|.down; 
oo, nd[q|.down =r, nd[r].up = q; 


This code is used in section 18. 


22. When we have input all options, we can tighten the upper bound of an item if there are less options 
available for that item than its upper bound. 


( Tighten bounds of the primary items 22) = 
for (k =1; k < second; k++) { 
if ((00, nd[k].len < cl[k].ubound)) { 
if (vbose & show_details ) 
fprintf (stderr, "Tightening, upperjboundforitem,"O".8syfrom,"O"dytoy"O"d\n", 
cl[k].name, cl[k].ubound, nd[k].len); 
o, cl[k].ubound = nd|k].len; 
bounds_tightened + ; 
} 
} 


This code is used in section 1. 


23. (Report the successful completion of the input phase 23) = 
forintf (stderr, "("O"11ld options, ,"O"d+"O"d items, .,"O"d entries,successfully read) \n", 
options, second — 1, last_itm — second, last_node — last_itm); 


This code is used in section 1. 
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24. The item lengths after input should agree with the item lengths after this program has finished— 
unless, of course, we’ve successfully simplified the input! I print them (on request), in order to provide some 
reassurance that the algorithm isn’t badly screwed up. 


(Report the item totals 24) = 
{ 
fprintf (stderr, "Item ,totals:"); 
for (k =1; k < last_itm; k++) { 
if (k = second) fprintf(stderr ,"|"); 
forintf (stderr, "\"O"d", nd[k].len); 


} 


fprintf (stderr, "\n"); 


} 


This code is used in section 1. 
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25. The dancing. The logic of preprocessing from DLX-PRE does not easily generalize to items with 
higher multiplicities. This program does not yet try to be very smart, and will often skip items with higher 
multiplities while preprocessing. A future version may potentially also try to handle those items. 

Let’s call a primary with lbound = ubound = 1 a primary item without multiplicities. Suppose p is a 
primary item with lbound > 0, and c is a primary item without multiplicities or a secondary item, such that 
every option containing p also contains an uncolored instance of c. Then we can delete item c, and every 
option that contains c but not p. For we’ll need to cover p,and then c will automatically be covered too. 

More generally, if p is a primary item with /bound > 0 and r is an option such that p ¢ r but every option 
containing p is incompatible with r, then we can eliminate option r: That option can’t be chosen without 
making p uncoverable. 

This program exploits those two ideas, by systematically looking at all options in the list for item c, as ¢ 
runs through all primary items without multiplicities and all secondary items. 

This algorithm takes “polynomial time,” but I don’t claim that it is fast. I want to get a straightforward 
algorithm in place before trying to make it more complicated. 

On the other hand, I’ve tried to use the most efficient and scalable methods that I could think of, consistent 
with that goal of relative simplicity. There’s no point in having a preprocessor unless it works fast enough 
to speed up the total time of preprocessing plus processing. 


26. The basic operation is “hiding an item.” This means causing all of the options in its list to be invisible 
from outside the item, except for the options that color this item; they are (temporarily) deleted from all 
other lists. 

As in DLX2, the neat part of this algorithm is the way the lists are maintained. No auxiliary tables 
are needed when hiding an item, or when unhiding it later. The nodes removed from doubly linked lists 
remember their former neighbors, because we do no garbage collection. 
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27. Hiding is much like DLX2’s “covering” operation, but it has a new twist: If the process of hiding item 
c causes at least one primary item p with lbound > 0 to become empty, we know that c can be eliminated 
(as mentioned above). Furthermore we know that we can delete every option that contains c but not p. 

Therefore the hide procedure puts the value of such p in a global variable, for use by the caller. That 
global variable is called ‘stack’ for historical reasons: My first implementation had an unnecessarily complex 
mechanism for dealing with several primary items that simultaneously become empty, so I used to put them 
onto a stack. 


(Subroutines 8) += 
void hide(int c) 
{ 
register int cc, 1, r, rr, nn, uu, dd, t, k=0; 
for (0, rr = nd[c].down; rr > lastitm; o, rr = nd[rr|.down) 
if (0, >nd[rrj.color) { 
for (nn =rr+1; nnZ#rr; ) { 
0, uu = nd[nn].up, dd = nd[nn].down; 
0, cc = nd[nn].itm; 
if (cc <0) { 
nn = uu; 
continue; 


oo, nd[uu].down = dd, nd[dd].up = uu; 

o,t = nd(cc].len — 1; 

o, nd|cc].len = t; 

if ((=0A ce < second A (0, cl[cc].lbound)) stack = cc; 
nnt++; 


} 


28. (Subroutines 8) += 
void unhide(int c) 
{ 
register int cc, 1, r, rr, nn, uu, dd, t; 
for (0, rr = nd[c].down; rr > last.itm; o, rr = nd[rr|.down) 
if (0, >nd[rr].color) { 
for (nn =rr+1; nn Arr; ) { 
0, uu = nd[nn].up, dd = nd[nn].down; 
0, cc = nd|[nn].itm; 
if (cc <0) { 
nn = uu; 
continue; 
‘ 
o,t = nd{[cc].len; 
oo, nd[uu].down = nd[dd].up = nn; 
o, nd[cc].len =t +1; 
nn; 
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29. Here then is the main loop for each round of preprocessing. 


(Reduce the problem 29) = 
for (cc = 1; cc < last_itm; cc++) 
if (0, (cc > second A nd{cc].len = 0) V (ce < second A nd{cc].len < cl[cc].lbound)) 
(Take note that cc has not enough options 41); 
for (rnd = 1; rnd < rounds; rnd++) { 
if (vbose & show_choices) fprintf (stderr, "Beginning, round,."O"d:\n", rnd); 
for (change =0,c = 1; c < last_itm; c++) 
if (0, nd[c].len A (c => second V (0, cl{c].lbound) A (cl[c].ubound = 1))) 
(Try to reduce options in item c’s list 31); 
if (~change) break; 
} 


This code is used in section 1. 


30. (Global variables 2) += 


int rnd; /* the current round */ 
int stack; /* a blocked item; or top of stack of options to delete */ 
int change; /*x have we removed anything on the current round? */ 


31. In order to avoid testing an option repeatedly, we usually try to remove it only when c is its first 
element as stored in memory. 

Note (after correcting a bug, 02 January 2023): If c is secondary and has a nonzero color in option r, we 
should not try to remove r, because r has not been hidden by the hide routine. Thus we might miss some 
potential deletions. Users can avoid this by putting all of the colored secondary items last in every option. 


(Try to reduce options in item c’s list 31) = 
{ 

if (sanity_checking) sanity(); 

if (delta A (mems > thresh)) { 
thresh += delta; 
forintf (stderr, "after "O"1ldmems:"O"d."O"d, "O"dyitemsout,"O"d options,out,", 

mems, rnd, c, itms_out, options_out ); 

forintf (stderr, "\,"O"d bounds tightened\n",, bounds_tightened ); 


if (mems > timeout) goto finish; 
stack = 0, hide(c); 
if (stack) (Remove item c, and maybe some options 32) 
else { 
for (o,r = nd[c].down; r > lastitm; o,r = nd[r].down) { 
for (¢=r-—1; 0,nd[q|.down =q—1; q—) ; /* bypass null spacers */ 
if (0, nd[g].itm < 0A (0, >nd[r].color)) 
/* ris the first (surviving, uncolored) node with multiplicity 1 in its option */ 
(Stack option r for deletion if it leaves some primary item uncoverable 36 ); 


unhide(c); 
for (r = stack; r; r=rr) { 
00, rr = nd[r].itm, nd[r].itm = c; 
(Actually delete option r 40); 
} 
} 
} 


This code is used in section 29. 
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32. (Remove item c, and maybe some options 32) = 
{ 

unhide(c); 

if (vbose & show_details ) 
fprintf (stderr, "Deleting item,"O".8s, forced by,"O".8s\n", cl[c].name, cl[stack].name); 

for (0,r = nd[c].down; r > lastitm; r= rrr) { 
0, rrr = nd[r].down; 
(Delete or shorten option r 33); 


o, nd|c].up = nd[c].down = c; 
0, nd[c].len = 0, itms_out +; /* now item c is gone */ 
change = 1; 

} 


This code is used in section 31. 


33. We're in the driver’s seat here: If option r includes stack, we keep it, but remove item c. Otherwise 
we delete it. 


(Delete or shorten option r 33) = 
{ 
for (q=r+1l,q#4r; ) { 
0, cc = nd[q].itm; 
if (cc <0) { 
0,q = ndlq].up; 
continue; 
} 
if (cc = stack) break; 
qtr; 
if (¢q#r) (Shorten and retain option r 34) 
else (Delete option r 35); 


} 


This code is used in section 32. 


34. (Shorten and retain option r 34) = 
{ 
if (vbose & show_details) { 
fprintf (stderr, "\jshortening,."); 
t = dpoption(r, stderr), fprintf (stderr, "\,Coption,,"O"d) \n",t); 
} 
0, nd[r].up =r +1,nd[r].down =r —1; /* make node r into a spacer */ 
0, nd[r].itm = 0; 


} 


This code is used in section 33. 
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35. (Delete option r 35) = 
{ 
if (vbose & show_details) { 
forintf (stderr, "Wdeleting,."); 
t = dpoption(r, stderr), fprintf (stderr, "\,Coption,,"O"d) \n",t); 


options.out ++; 
for (o,g=r+1; qr; ) { 
0, cc = nd{qj.itm; 
if (cc <0) { 
0,q = nd{q].up; 
continue; 


o,t = nd(cc].len — 1; 
if ((cc > second A at) V (cc < second A (0,t < cl[cc].lbound))) 
(Take note that cc has not enough options 41); 
o, nd[cc].len = t; 
o, uu = nd{q|.up, dd = nd[q|.down; 
oo, nd[uu|.down = dd, nd[dd].up = uu; 
qty; 
} 
} 


This code is used in section 33. 
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36. At this point we’ve hidden item c and option r. Now we'll hide also the other items (if primary, only 
when ubound = 1) in that option; and we’ll delete r if this leaves some other primary item uncoverable. (As 
soon as such an item is encountered, we put it in pp and immediately back up.) 

But before doing that test, we stamp the aux field of every non-c item of r with the number r. Then we’ll 
know for sure whether or not we’ve blocked an item not in r. 

When cc is an item in option r, with color x, the notion of “hiding item cc” means, more precisely, that 
we hide every option in cc’s item list that clashes with option r. Option rr clashes with r if and only if 
either x = 0 or rr has cc with a color ¥ x. 


(Stack option r for deletion if it leaves some primary item uncoverable 36) = 
{ 
for (g=r+1;;) { 
0, cc = nd{q].itm; 
if (cc <0) { 
0,q = nd[q).up; 
if (¢>7r) continue; 
break; /* done with option */ 
} 
o, nd[cc].aux = r,q+t+; 
} 
for (pp =0,g=rt+1;; ) { 
0, cc = nd[q].itm; 
if (cc <0) { 
0, q = nd{qj.up; 
if (¢>1r) continue; 
break; /* done with option */ 
} 
if (cc > second V (o, cl[cc].ubound = 1)) { 
for (x = nd[q].color,o0,p = nd[cc].down; p > last_itm; o,p = nd[p|.down) { 
if (x > 0A (0, nd[p].color = x)) continue; 
(Hide the entries of option p, or goto backup 37); 


qt; 


backup: for (q=r—1; q4r; ) { 
0, cc = nd{qj.itm; 
if (cc <0) { 
0,q = nd[q|.down; 
continue; 
} 
if (cc > second V (o, cl[cc].ubound = 1)) { 
for (x = nd[q].color,o,p = nd[cc].up; p > lastitm; 0,p = nd[p].up) { 
if (c > 0A (0, nd[p].color = x)) continue; 
(Unhide the entries of option p 38); 
} 
} 
q—; 
} 


if (pp) (Mark the unnecessary option r 39); 


} 


This code is used in section 31. 
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37. Long ago, in my paper “Structured programming with go to statements” [Computing Surveys 6 
(December 1974), 261-301], I explained why it’s sometimes legitimate to jump out of one loop into the midst 
of another. Now, after many years, I’m still jumping. 
(Hide the entries of option p, or goto backup 37) = 
for (qg=p+1; q@#p; ) { 
0, cc = nd[qq].itm; 
if (cc <0) { 
0, qq = nd[qq].up; 
continue; 


o,t = nd[cc].len — 1; 

if (cc < second A nd{cc].aux #r A (0,t < cl[cc].lbound)) { 
pp = ce; 
goto midst; /* with fingers crossed »*/ 

} 

o, nd{cc].len = t; 

o, uu = nd{qq].up, dd = nd[qq].down; 

oo, nd[uu].down = dd, nd[dd].up = uu; 

qq; 


} 


This code is used in section 36. 


38. (Unhide the entries of option p 38) = 
for (qq=p-1; q@#p; ) { 
0, cc = nd[qq].itm; 
if (cc <0) { 
0, gq = nd[qq|.down; 
continue; 


00, nd|cc].len +; 

o, uu = nd{qq].up, dd = nd[qq].down; 

oo, nd[uu].down = nd[dd].up = qq; 
midst: qq--; 


} 


This code is used in section 36. 
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39. When I first wrote this program, I reasoned as follows: “Option r has been hidden. So if we remove 
it from list c, the operation unhide(c) will keep it hidden. (And that’s precisely what we want.)” 

Boy, was I wrong! This change to list c fouled up the unhide routine, because things were not properly 
restored/undone after the list no longer told us to undo them. (Undeleted options are mixed with deleted 
ones.) 

The remedy is to mark the option, for deletion later. The marked options are linked together via their 
itm fields, which will no longer be needed for their former purpose. 


(Mark the unnecessary option r 39) = 


if (vbose & show_details) { 
forintf (stderr, "\"O" .8syblocked by", cl[pp].name); 
t = dpoption(r, stderr), fprintf (stderr ,"\,Coption,"O"d) \n", t); 


options.out++, change = 1; 
0, nd[r].itm = stack, stack =r; 


This code is used in section 36. 


40. (Actually delete option r 40) = 
for (p=r+1;; ) { 
0, cc = nd|p].itm; 
if (cc <0) { 
0, p = nd[p].up; 
continue; 
} 
0, uu = nd[p].up, dd = nd[p|.down; 
oo, nd[uu].down = dd, nd[dd].up = uu; 
oo, nd[cc].len—-; 
if ((cc > second A nd{cc].len = 0) V (cc < second A (0, nd[cc].len < cl{cc].lbound))) 
(Take note that cc has not enough options 41); 
if (p =r) break; 
pr; 
i 


This code is used in section 31. 


41. (Take note that cc has not enough options 41) = 
a 
itms_out ++; 
if (cc > second) { 
if (vbose & show_details) fprintf (stderr ,"4"O" .8s,has, mot enough, options, ("O"d<"O"d) \n", 
cl[cc].name, nd[cc].len, cl[cc].lbound); 
} else (Terminate with unfeasible item cc 42); 


} 


This code is used in sections 29, 35, and 40. 
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42. We might find a primary item that has not enough options left. In such a case all of the options can 
be deleted, and all of the other items! 
(Terminate with unfeasible item cc 42) = 


if (vbose & show_details ) 
forintf (stderr, "Primary,item,"O".8s,has,notenough, options, ("O"dy< "O"d) !\n", 


cl[cc].name, nd[cc].len, cl[cc].lbound); 
options_out = options; 
itms_out = lastitm — 1; 
printf (""O".8s\n", cl[cc].name); /* this is the only line of output «/ 
goto alldone; 


} 


This code is used in section 41. 


22 THE OUTPUT PHASE DLX3-PRE §43 


43. The output phase. Okay, we’re done! 


(Output the reduced problem 43) = 
(Output the item names 44); 
(Output the options 45); 


This code is used in section 1. 


44. In order to be tidy, we don’t output a vertical line when all the secondary items have been removed. 


(Output the item names 44) = 
for (c=p=1; c< last_itm; cH+) { 
if (c= second) p=0; /* no longer primary */ 
if (0, nd[c].len) { 
if (p =0) p= —1, printf ("ul"); 
if (p > 0A ((o, cl[c].lbound F cl{c].ubound) V cl[c].ubound £1)) { 
if (cl[c].lbound & cl[c].ubound) 
printf ("\"O"d:"O"dl|"O".8s", cl[c].lbound, cl{c].ubound, cl[c].name); 
else printf ("y"O"d|"O".8s", cl[c].ubound, cl[c].name); 


else printf (""O".8s", cl[c].name); 


} 


i 
printf ("\n"); 


This code is used in section 43. 


45. (Output the options 45) = 
for (c= 1; c< last_itm; c++) 
if (0, nd[c].len) { 
for (0,r = nd[c].down; r > lastitm; o,r = nd[r].down) { 
for (¢=r-—1; o,nd[q|.down =q—1; q—) ; 
if (0, nd[q].itm <0) { /* r was the leftmost survivor in its option */ 
t = dpoption(r, stdout); 
printf ("\n"); 
if (vbose & show_orig_nos) printf ("| (from ,"O"d) \n", t); 
} 
} 
} 


This code is used in section 43. 
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